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We consider weak solutions to dispersive partial differential equations with periodic boundary con- 
ditions and initial data with jump discontinuities. These are already known to be continuous at 
irrational times and piecewise constant at rational times; we show that as time approaches a ratio- 
nal value the solution exhibits a ringing effect, with the characteristic overshoot of fixed amplitude 
near the discontinuities. Furthermore this effect is the same whether the sequence of times follows 
Qh ' rational or irrational values. 
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1. Introduction 

Consider the equation 

U^ n \t,x) = {2%if- l U { ^\t,x) (1.1) 

for n>2. Any classical solution to this equation satisfies the condition that 



ffu(t,x)^°' n \t,x)-(-2ni) n - 1 F^' \t,x)^dxdt = (1.2) 



for any ff°° function F of compact support in M 2 , as can be seen using integration by parts, 
and we call a function U(t,x) satisfying (1.2) a weak solution to (1.1). We are interested 
here in weak solutions with periodic boundary conditions U(t ,x+ 1) = U(t ,x) for all x. It is 
natural to study such solutions through their representations by Fourier series, which are of 
course ubiquitous in mathematics, and we will use techniques from analytic number theory 
to analyse these series and show some asymptotic properties of the weak solution U(t,x). 

To motivate this, for the present consider the simpler case of (1. 1) for n = 2 with vanish- 
ing boundary conditions U(t,x) — > as x — > ±°°, and initial data with jump discontinuities; 
for instance U(0,x) = %{x), where % is the characteristic function of the interval [— y,y] for 
< y < 1/2. By considering the Fourier transform in x of U(t,x) one can show that 

sm2nyk 



U(t,x)= J U(0,k)e(tk 2 +xk)dk= f & -^-y-e{tk 2 +xk)dk (1.3) 
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where we have used the notational convention e(t,) = e 2m ^, as throughout this paper. With 
t = 0, the integral represents the initial data %(x) in L 2 . Moreover it is not hard to prove that 
the truncated integral 

rN sm2nyk 



-^e(xk)dk 

-N 7lK 



converges pointwise to %(x). Of course, the convergence cannot be uniform and indeed one 

has 

1 1 f&mlnw 

lim U (0,y+ sN;N ) = --- / dw. 

n^oo 2 n Jo w 

This is the well-known Gibbs phenomenon, which is an oscillatory ringing produced by 
truncation of the Fourier representation of a function with a discontinuity. 

There is, however, an entirely different phenomenon, which is an oscillatory ringing 
produced not by truncation but rather by a different mechanism which, as we will explain, 
should be thought of as a dispersive regularization of a discontinuity. Indeed, by contour 
deformation or otherwise, the integral solution (1.3) can be seen to be analytic in t or x if 
t > 0. Away from x = ±y the integral converges nicely to U(0,x) as t — > + , however in 
the rescaled variable s where x = zty+st 1 ' 2 one can obtain the asymptotic expression (see 
DiFranco & McLaughlin 2005) 

U (t,y+ S t^) = \ - ^Erf (^1 *"*/*,) + O (t" 2 ) (1.4) 

where Erf denotes the usual error function (see Abramovitz & Stegun 1972). Similar phe- 
nomena can be shown to occur for n > 2, although with different ringing functions. 

Our goal here is to study similar effects in the more complicated situation where U(t,x) 
obeys periodic boundary conditions, which have not been observed previously, and to give 
asymptotic expressions for these effects. Here Fourier series expansions replace the Fourier 
transform (or if one prefers, series replace integrals), and the resulting structure is more 
complicated, presenting ringing effects at all rational times, not just at zero as in the "whole 
line" case. The general setup is as follows. Suppose now that U(t,x) is defined on [0, 5]xK 
for some 8 > and is periodic in x, that is, U(t,x) = U(t,x+ 1) for all t and x. Such a 
function has a Fourier series 

l 
U(t,x) ~ £ c k {t)e(kx) , c k (t) = fu(t&)e{-kQa%. (1.5) 

*=— 

If U(t,x) is a square-integrable weak solution to (1.1) then standard arguments show that 
Ck(t) = Cke{tk n ) where the constants q are the Fourier coefficients of U(0,x), and where 
the convergence of the series to U(t,x) is understood in the L 2 sense. Note that periodicity 
in x has forced periodicity in t also, which is the root cause of the special behaviour at 
rational times. We will say a function f{x) is in class & if it is integrable, periodic of 
period 1, piecewise continuously differentiable, and f{x) = {f(x + ) +f(x~)}/2 for all x. 
It is a well-known theorem of harmonic analysis (see Katznelson 2004, for instance) that 
functions / of class @ have Fourier series that converge pointwise to /, in the sense that 

f(x) = lim £ c k e(kx) , 

°°\k\<K 
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that is, the series converges pointwise once written in terms of sines and cosines rather than 
exponentials. (The series of exponentials may diverge, in the literal sense). We will consider 
periodic discontinuous initial conditions 

U(0,x)= £ X[-y,t](* + (1-6) 

where < y < 1/2, which is in class % with q = (71&) _1 sin27tyfc, so the Fourier series 
(1.5) for 17 (f,*) is 

C/( f , x ) „ V *^t e (tk n +xk) . (1.7) 

The series clearly converges in the L 2 sense, and represents the unique weak solution in that 
space. We will show some structural properties of U(t,x) by identifying it with this series 
representation and considering the properties of the series. 

To understand why the situation should be so much more complicated requires review- 
ing some known results. In a variety of applied sciences, the term "dispersive" describes 
a localized quantity which spreads out as time passes. In the context of wave phenomena, 
an equation is called dispersive if its simple oscillatory solutions propagate at velocities 
that depend on their spatial frequencies. This is the case here; rewriting the coefficients as 
e{k{x + k n t)) they can be seen as a family of traveling waves with velocity k"~ , which 
is clearly dependent on the spatial frequency k. Thus if one starts with localised initial data 
the various terms of the Fourier series of the solution spread out in space, since they move at 
different velocities. As time passes the solution will appear more extended, with the highest 
spatial frequency components most apparent in the furthest reaches of the solution. Prop- 
erly speaking, dispersion requires sufficient extent to permit the observation of spreading. 
This may, or may not, be the case for a partial differential equation with periodic boundary 
conditions, but nonetheless such an equation is still called dispersive if its family of solu- 
tions have frequency dependent velocity. In this context it is striking that periodic boundary 
conditions should cause a dispersive equation like (1.1) to show recurrence, as well as rad- 
ically different behaviour at rational and irrational times. This can be seen experimentally 
in the case n = 2, as reported by Talbot 1836, who shone light on a periodic grating and 
looked at the images produced by it. He reports "...a regular alternation of numerous lines 
or bands of red and green colour, having their direction parallel to the lines of the grating. 
On removing the lens a little further from the grating, the bands gradually changed their 
colours, and became alternately blue and yellow." (Talbot 1836). This is now known as the 
Talbot effect, and has been extensively studied; we refer the reader to Berry & Klein 1996 
for details of the effect, but also for two contributions which are pertinent to our discussion 
here. Supposing that the period of the grating is a and that the wavelength of the light is 
X, the images seen by Talbot form at regular multiples of a 2 fk. (In the context of (1.1), 
we can view Talbot's grating as initial data which is periodic, being 1 on the slits of the 
grating and zero elsewhere). Berry and Klein observe that at rational multiples (p/q)cr/X 
of this distance " ...these fractional Talbot images consist of q equally spaced copies of the 
transmission function of the grating, which superpose coherently when they add up.". Fur- 
thermore, they show that these translates have phases given by Gauss sums, so the solutions 
can be written down explicitly at rational times, and are piecewise constant functions of x. 
(This was also discussed in Olver 2010). They also observe that this is in sharp contrast 
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with irrational times; considering the solution at a sequence of rational times tending to an 
irrational t, they show that "the graph of a function with power spectrum \g n \ 2 proportional 
to n~v is a fractal curve with fractal dimension D = (5 — P)/2. A smooth curve has D = 1, 
and a curve with D = 1 is so jagged that it is almost area filling; curves with 1 < D < 2 
are continuous but non-differentiable. Thus, since thefourier coefficients of the initial data 
have power spectrum decaying like n~ 2 , the fractal dimension must be 3/2." 

A number of researchers from a variety of different areas of analysis, including Arkhipov 
& Oskolkov 1989, Stein & Wainger 1990, Oskolkov 1992, Kapitanski & Rodnianski 1999, 
Rodnianski 1999 and Rodnianski 2000, put these ideas on a rigorous footing by establishing 
convergence properties of the associated Fourier series in the space of continuous functions. 
One result which is particularly pertinent to the discussion here is due to Rodnianski 2000: 
this discusses periodic solutions to (1.1) with n = 2 with initial data of bounded variation 
but not in the space U £> oH l ' 2+£ . (Here the Sobolev space H s is the space of functions 
whose Fourier coefficients f(k) satisfy YX=\ l/WI 2 (l + |^| 2 ) < °°.) The result states 
that for almost all irrational t, including all algebraic t, the fractal dimension of the graphs 
of the real and imaginary parts of U(t,x) (which are continuous functions of x) is 3/2. 

In fact there is a dichotomy between rational and irrational times for all n > 2, which 
is reflected in our calculations. In contrast to the considerations in Berry & Klein 1996 
or Rodnianski 2000, which consider sequences of rational times tending to an irrational 
limit, to understand the nature of the function at irrational times, we consider sequences of 
times, rational or irrational, tending to a rational limit, to understand the ringing phenomena 
exhibited. 

As observed above (1 .7) converges in L 2 to U(t ,x), however convergence in any stronger 
sense is more complicated, depending greatly for the reasons mentioned above on whether t 
is rational or irrational. For irrational times t we will use rational approximations; we speak 
of approximants as being reduced fractions u/q such that 



u 
t 



< 3 • d" 8 ) 



There are infinitely many approximants to any irrational t, some of them given by the con- 
vergents Uj/qj € Q from the continued fraction expansion to t, which obey the inequality 



u, 
t - 



4j 



< — *— ■ (1.9) 

qjqj+i 



(See Hua 1982 or Hardy & Wright 1979). Note also that a celebrated theorem in Roth 
1955 states that if t is an algebraic irrational then for any r\ > there are only finitely many 
solutions to 

We need to introduce some sets of irrationals described by their rational approximations. 
To this end n > 2 will always represent the order of the equation (1.1), and A will be fixed 
once and for all in (0, 1), and in (0, 1/2) in the case n = 2. Note that several of the sets and 
implied constants below will depend on n and A without this being explicitly mentioned or 
shown in the notation. 
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Definition 1.1. Define srfto be the set of t G [0, 1] such that for all sufficiently large M 
there is an approximant u/q to t with 

M A <q<M n ~ A (1.11) 

and define si m to he the set oft £ [0, 1] such that for all M >m there is an approximant u/q 
to t satisfying (1.11). Further, for OC >0 define S3 m & to he the set of t € [0,1] suchthatfor 
all M > mt~ a there is an approximant u/q to t satisfying (1.11). 

Although some parts are not strictly necessary below, we prove the following lemma. 

Lemma 1.2. Let /u denote Lebesgue measure, and suppose (n — A) -1 < a < A . 

1. If mi < m2 then srf m C srf mi . 

2. Ift € stfthen there exists m such that t G stf m , so srf= \J m >2- s ^m- 

3. There is a constant c n ^ such that /u (£/ m ) > 1 — c n A>n i+2A ~ n , which tends to 1 as m — > oo, 
so n{sd) = 1. 

4. ^contains all algebraic irrationals. 

5. For all m, stf m C £% m ^, so fi (S2 mi a) —■ 1 as m —> oo. 

6. There is a constant c n ^ a such that 

,(«-l-2A)a m _„ +A +l 



A< {^m,a n [0,f ]) > to 1 - C„ Aa t 



The convergence of the series (1.7) is described in the following result. 
Theorem 1.3. Let 

Sk= e :<*£*>. (U2) 

o<|-t|<^ 

1. Ift is rational then Sk converges pointwise in x to a piecewise constant function. 

2. Ift £ s^then the sequence Sk converges uniformly in x. 

Note that Arkhipov and Oskolkov 1989 have shown that a more general class of series 
converges pointwise, and hence the partial sums (1.12) converge pointwise in x for all 
times. Their proof uses Vinogradov's method for exponential sums in place of the Weyl 
shift method described below and used here; we state and prove Theorem 1.3 since it is 
necessary to Theorem 1 .5 below, for which we are as yet unable to apply the Vinogradov 
method. 



As discussed above, Part 1 of this Theorem has been noted elsewhere. The same result 
is true, however, in a much more general context than we can prove Part 2, so we state and 
prove it separately, and deduce 1 of Theorem 1 as an immediate corollary. 

Theorem 1.4. Let P be a polynomial with integer coefficients and L be the differential 
operator 

L = 2nip(±l 

\ 1%i ox 

and consider the initial value problem 

LU(t,x) = U t (t,x) (1.13) 
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Figure 1: The real part of U(t,x) torn = 2 at u/1 for u = 0,1,2,3 and u = 4,5,6,7 
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Figure 2: The imaginary part of U(t,x) for n = 2 at u/1 for u = 0, 1,2,3 and u = 4,5,6,7 



with U(0,x) = f{x), where f G & and f{x+ 1) = /(x). Denoting 

G(u,v;q)= £ e q {uP{w) - vw), 

w(mod <?) 

a? ?/je rational time t = u/q we have 

U(t,x) = - £ G(«,v;<?)/(x+-J . 

r/jwi ifU(0,x) is piecewise constant, U(t,x) is piecewise constant at all rational times. 

We can rephrase this theorem as saying that at rational times t the solution to the initial 
value problem ( 1 . 13) is a linear combination of translates of the initial data. Note also that 
very generally (see Schmidt 2004) we have G <C q 1 ' 2+E , a bound which is not necessary to 
prove Theorem 2, but which we will use elsewhere; in particular note that 



(1.14) 



(1.15) 



£ e q {uw n -vw)<ZL{u,v,q) l l 2 q l l 2 ^ 

w(modg) 



(1.16) 
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where (u,v,q) denotes the greatest common divisor of u,v and q, and where the implied 
constant depends on £ and n. As noted in Berry & Klein 1996, in the special case of the 
Schrodinger equation these sums are Gauss sums and can be evaluated for general q, but in 
general this is not realistic and we must be satisfied with estimates. We omit the proof of 
(1.16), merely noting that it follows from Theorem 2.5 in Schmidt 2004, which contains an 
extensive discussion of these issues. 




mk0Aj \J^)0<m 
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Figure 3: U(t,x) for n = 2 and u/q = 1/17,1/33,1/65, 1/129.1/257,1/513, and 
1/2049,1/4097,1/7374 



An example is given by the solution to (1.1) in the case n = 2 and y= 1/tc, at the ratio- 
nal times t = u/1. Its real and imaginary parts are plotted in Figure 1 and 2, and can be seen 
to start and finish equal to the initial data. Note that for odd n, of course, the solutions are 
purely real. Note also that Figures 1-6 and 10 were produced with the aid of the explicit 
closed form expressions afforded by Theorem 1.4, in all cases for rational values of t. In- 
deed, formula (1.15) represents the solution U(x,t) as a finite sum. The ringing effect which 
is the main topic of this paper occurs as t tends to a rational value. To illustrate this effect 
in Figure 3 we plot the real part of U(t,x) at rational times \/q for various denominators 
q. The first values u/q = 1/17, 1/33, 1/65 are too large for any asymptotic tendency to be 
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Figure 4: Zoomed-in sections from U(t,x) for n = 2 and u/q = 1/1374 



evident, although they clearly show the piecewise constant behaviour described in Theorem 
1.4, and the number of segments grows approximately linearly with q. The initial data starts 
to appear for u/q = 1/129, 1/257, 1/513, and the characteristic overshoot of a ringing ef- 
fect becomes evident for u/q = 1/2049,1/4097,1/7374. Of course, in these last figures 
the resolution is insufficient to show the piecewise constant behaviour, which is clear in 
zoomed-in plots from the case q = 131 A shown in Figure 4. The ringing is also clear in the 
case of odd n, which shows some features similar to that of even n and others which are 
quite distinct. In particular, the solutions are not even as functions of x, although they are 
purely real, and the overshoot is asymmetrical on the two sides of the discontinuities. This 
can be seen in Figure 5, which graphs the solution for n = 3, again with y = 1/71 at several 
times \/q. 

Figures 1-6 and 10 were produced using Matlab. Fixing n and q, we first produce a 
q x q matrix containing the values of G(u,v;q). Then a vector representing the initial data 
f(x) is produced for a fixed fixed uniform grid of points x between —1/2 and 1/2, and the 
sum in (1.15) evaluated as a matrix- vector multiplication. Matlab calculates using IEEE 754 
double-precision format (1 sign bit, 1 1 exponent bits, 52 mantissa bits), which translates into 
approximately 16 decimal digits of accuracy. For our purposes, this is sufficiently accurate; 
in fact the real limitation is not the precision but rather the time required to compute the 
matrix G, and we are limited by this restriction to roughly q f=a 100000. 

The main result of this paper is the following theorem, which confirms the two principal 
aspects of these plots: that the pointwise limit is equal to the initial condition, and that there 
is a ringing effect, given by an overshoot of fixed amplitude which migrates towards the 
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discontinuity. 

Theorem 1.5. LetU(t,x) be the solution to (1,1) with periodic initial conditions U(0,x) 
given by (1.4). Taking any m>2 and any OC between l/(n — A) and l/(n — 1), let t — > + 
following a sequence of points such that t = u/q with t < mq '°^ n ' if rational, and 
t 6 £§m t a if irrational. Then for any £ > 0, 

l + (9^aA/2"- 1 - E + ^| x 2_ y 2|-A ,j| x | <y 
U(t,x) = < / V , N J 

0[t aA / 2 - £ + t a \x 2 -f\- 1 ) if\x\>J 



and ifS > is fixed then for s £ [—S,S], 

oo 

U (t,± 1+ st l l n ) =\t^ J 'e(y n )sm2Ksy^ + 0(t aA / 2 "- l -*+t m - 1 

— oo 

for even n and 

oo 

u(t,±y+st l ^=^T^fe(y n )sm2n(y n + sy)^ + 0(t aA ^ l - £ + t 



arc-1 



for odd n, where the implied constants depends on S, J, m and £. 

Note that it is the dependence of the implied constant on m that forces us to consider 
t tending to zero in the fixed set <^ m , a rather than in the union of all such, which is stf. 
Nonetheless, there are plenty of points in the set £$ m>a , as shown in part 6 of Lemma 1 .2. The 
restriction on rational values of q is essentially a requirement that the numerator in t = u/q 
not remain large; that is, that as q grows, t tends to zero reasonably quickly. Note that it 
is similar to the definition of the set -SS mM in the irrational case, although less restrictive 
on q. In fact, in the rational case the bounds proved are a little stronger than stated in the 
Theorem, as can be seen from the proof. 

The ringing effect is described in Theorem 1.5 in a renormalised variable s; it is in- 
teresting to graph this renormalised behaviour; in the case n = 2 Figure 6 shows the real 
and imaginary parts (left and right, respectively) of U(t,s + yt 1 ^ 2 ) for t = 1/202, 1/1616, 
1/6464, and 1/51712. By way of comparison, we can also plot the integral to which U 
is asymptotic; for n = 2 the real and imaginary parts are shown in Figure 7, which appear 
very close to the plots of U for large denominator. Although it is computationally difficult 
to produce plots of the function U{u/q,x) with large q for n > 3 larger than 3, we can plot 
the integral from Theorem 1.5, which is much easier to calculate. (This is, of course, the 
whole point of proving asymptotic expressions!) For instance, Figure 8 shows the real and 
imaginary parts of the integral for n = 6. In the case of odd n we can produce similar plots, 
although as above they do not show the same symmetry as for even n; for instance we have 
the graphs for n = 3 and n = 17; Figures 7-9 were produced by Mathematica. In fact the in- 
tegrals appearing in Theorem 1.5 are, for each integer n>2, special functions representable 
using the hypergeometric function p F q (typically with p = 1 or p = 2). 

As mentioned above, while Theorem 1.5 describes the ringing effect in detail for times 
tending to zero, in fact the same phenomenon is repeated at each of the jump discontinuities 
for rational time u/q, as t — > u/q. This follows immediately from Theorems 1.4 and 1.5 
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considered together. An example is shown in Figure 10, which shows the solution for n = 2 
and y = 1/51. at t = 468/3277, which is a close approximation to t = 1/7. Note that the 
graph exhibits ringing effects near each of the jump discontinuities in the graph of U for 
t = 1/7 as in Figure 1. 
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Figure 5: « = 3 and u/q = 1/65, 1/257, 1/1025 and 1/65, 1/257, 1/1025 



As mentioned above the proof of Theorem 1.5 uses well-known estimates for Weyl 
sums which are obtained in the literature using the Weyl shift method; we state the required 

results before proceeding. If f(k) = a n k n -\ \- do, (X n ^ 0, is a polynomial of degree n > 2 

with real coefficients, and k varies over an interval of at most /n consecutive integers, then 
for any £ > 

\^e(f(k))\ N «^+ fJ N - £ min {^ ! 

l<n,-,>v>-i<ju-l l> 



«^'-'+/-' H 



AT-l 



l<m<G"-l) 



2{n\a n r\r 2 ...r n -i} 
1 



V min < u, , 

*- „_, Y'2{n\a n m} 



(1.17) 



where {x} denotes the distance from x to the nearest integer, N = 2 n l . (Originally due to 
Weyl. See Titchmarsh 1986). Note also that with a„ = u/q we have 

Y^ min^/7, — — ->= Y, rmn- li u i^ _ f Y 1 

l<m<(u-l)»- ] ^ 2{n.a n m\ J v{modq) < 2v i< m « u -i)"-' 



l<m<0u-l)" 

n\um=v(modq) 



<(//" l q l + l)(M+qlogq) 



and hence 



\£e(f(k))\ N «^ ! +v N -" +£ (v n - l q- 1 + l)G«+«)log ? • 



(1.18) 
(1.19) 



Furthermore, if 



h 
a n -~ 

<7 



< — > (&>?) = 1 > 
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Figure 6: Real and imaginary parts of U(t,y + st 1 ' 2 ) for f 

1/202, 1/1616, 1/6464, 1/51712. 
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Figure 7: The real and imaginary parts of the integral for n = 2. 




Figure 8: The real and imaginary parts of the integral for n = 6. 




Figure 9: The integral for n = 3 and n = 17. 



then for any £ > 0, 



M 

I 

m= 1 



£ e(f(m)) « {M- l +q- l +qM-"} M l+e q s 



(1.20) 



where the implied constant depends on 8, and is uniform in the coefficients 0C„_i,. . . ,(Xo. 
(See Hua 1965). 



The structure of the paper is as follows. In section 2 we establish the properties of the 
sets £/, £/ m and 33 mM stated in Lemma 1.2. Some of these are not strictly necessary for 
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Figure 10: The graph of U (468/3277, x) 

the remainder of the paper, but item 6 in particular is important, since it establishes that 
there are plenty of points in 33 mM through which t may tend to zero. In section 3 we prove 
Theorem 1.4, which handles the convergence of the Fourier series at rational times, and in 
section 4 we use (1.20) to consider convergence at irrational times in jz/and prove Theorem 
1.3. The proof of Theorem 1.5 follows in sections 5 to 7. As is common with periodic 
problems, one would like to use harmonic analysis such as the Poisson summation formula 
to replace the series by an integral, which is more easily analysed. In this case, however, it is 
necessary to first separate into large and small frequencies k, and treat the large frequencies 
separately, using quite different methods for rational times and for irrational times in SS m ^- 
These sections also use (1.19) and (1.20), although the arguments are more involved than in 
the proofs of convergence, and require that irrational times be in the set 3§, lha rather than si. 
In section 7 we apply the Poisson summation formula to the small frequencies, and analyse 
the resulting integral, in essentially the same way for rational and irrational times, which 
completes the proof. 

2. Proof of Lemma 1.2. 

Items 1 and 2 are clear from the definition of the sets. To prove item 3, consider t ^ g/ m ; 
thus there exists M > m for which there is no approximant u/q to t with M A < q < M n ~ A , 
and hence there must be consecutive convergents Uj/qj and Uj+i/qj+\ for which qj < M A 
> M n ~ A for some j. Defining 



and qj 



+ i 



ft 



q,M 



0, 



1 



qM n ~ 



uU 



1 



u 

- + ■ 



1 



q qM n - A q qM n 



U 



1 



1 



for any q and M, by (1.9) we have 
1 



t G 



0, 



qjqj+i 



u=\ 



1 



1 



- + ■ 

qj qj q j+l qj qjq j+ i 



U 



qM n - 



1 
qjqj+i 



,1 



and hence t € y q ,M for some q,M with q < M A and M >m. Thus 

(e[J U ^iM 

M>m q <M A 

and since /i (<5^,m) = M~ n+A , we can bound the measure of this union by 



-» i u u ■*** 

M>m q <M A 



< L M 

M>m 



2A -"<m 1+2A -". 
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This proves item 3. 

Suppose now that t is an algebraic irrational, and recall (1.8) and (1.10). Since all but 
finitely many convergents Uj/qj differ from t by at least q- T| , but by at most qjqj +l , it 

follows that for all but finitely many convergents we have qj + \ < q • . Now any positive 
number greater than 1 must fall between some qj and qj+\, hence for any M > 1 there 
will exist j for which q } < M n ~ A < q j+l < q)^ and hence M {n -^K l+ ^ < qj < M"~ A . It 
follows that for all sufficiently large M we can find a fraction uj/qj satisfying the conditions 
of Theorem 1 as long as we choose r\ < n/A — 2. This proves item 4. 

Item 5 is again clear from the definitions; item 6 can be proved by adapting the argument 
for item 3; if < t < to and t ^ 3§ m ,a there must exist M > mt~ a such that there is no 
approximant u/q to t with M A < q < M n ~ A . Since there is certainly an approximant with 

1 
< qQ 

where q < Q and Q is the integer part of M"~ A , it must be that q < M A . Were u = 
we would have t < {qQY 1 , which implies t < 2/(mr a ) n ~ A since Q > M n ~ A /2. Thus 
t < 2m- n+A t< n - A ^ and hence r«("- A )- 1 > m"~ A /2, which is clearly false, so u ^ 0. Thus t 

is in the set 

u 2 u 2 



=^,mC (J 

\<U<Uq 



q qM n - A ' q qM"- A 



for some M > mt~ a > mt^ a and q < M A , where uq is such that to < uo/q + 2/qM"~ A . The 
smallest element of .^m is at least \/2q, so t > \/2q and hence 2qt$ > 1. Thus we may 
choose uq to be the integer part of 2qt$ plus one, which is no more than Aqt®. The measure 
of ^ q ,M is thus no more than \6to/M n ~ A , and the measure of the union is bounded by 

/<( U U^.m)< I L^«^oKT" +2A+1 . 

\M>m^ a q<M" j M>mt tt q<M A 

This proves item 6, and completes the proof of Lemma 1. 

3. Convergence at rational times and proof of Theorem 1.4 

This follows from Theorem 2, which we prove. Fourier modes of the form e(tP(k)+kx) 
solve the equation LU = U t , so given / £ @, if c/c are the Fourier coefficients of the initial 
condition U(0,x) = f(x) then the Fourier series 

Y J C k e{tP(k)+kx) 
k 

is a solution to the initial value problem in the L 2 sense (as discussed in the introduction), 
and we wish to show it is convergent at rational times. From Katznelson 2004, for instance, 
we know the Fourier series for f(x) is convergent, so 

lim - V G(u,v;q) Y c k e[ (x + - )k) 
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exists and is equal to the right-hand side of (1.15). On the other hand, opening G and 
changing orders of the two sums this is equal to 

lim - Y c k Y e( -P(w)-vw + (x+- )k) = lim V c k e [ -P(k) +xk 

which hence exists, and proves both the convergence and Theorem 1.4. 

4. Convergence at irrational times and completion of the proof of Theorem 1.3 

Suppose now that P(k) is a monic polynomial of degree n > 2 with real coefficients. If 
t e j/then there exists m such that t € jrf m , and we consider K > m and L> K + m, and sum 
by parts to obtain 

L 

I < tP <V +xk K \ £ e(tP(k)+xk) + f £ e{tP{k)+xk)% (4.1) 

K<k<L K L K<k<L J K K<k<t, ^ 

The interval K < £, < K + m contributes at most 0(K~ 2 ) to the integral (where the implied 
constant depends on m, and hence on t, although this is unimportant in this section). In the 
remaining range ^>I + mwe rewrite the sum as so that 

S:= £ e(tP(k)+xk)= £ e(P*(k)) (4.2) 

K<k<% k<%-K 

where shifting k has not affected the leading coefficient of P, so P* (k) is a polynomial of 
the form required for (1.19) and (1.20), with a n = t. Choosing an approximant Uj/qj to t 
with (§ - K) A < qj < (§ - iT)"- A we apply ( 1 .20) and find 

5<<{^-^)- 1 +^7 1 +^^-^)^} 2l '^-^) 1+e ^<<(^-^) 1 " 2l " A+( ' I+1)e - (4-3) 
Specifying e = A/2 n (n + 1) this gives S <C (£, - iT) 1-8 , where 8 = A/2" and hence 

L 



^ e[tk"+xk) r _ 5 /Vi-8 Jtt 

2* — i — < L + / 5 ^ 



K<k<L 



where we note that the constants are independent of x. It follows that the sequence of partial 
sums (1.12) is a Cauchy sequence, and the series is thus convergent. The uniformity follows 
since none of the constants appearing depend on x, and hence the tail of the series for k > K 
is uniformly bounded by 0(K ). 

5. Bounding the contribution from large frequencies at irrational times 

Let t € & m ,a as i n Lemma 1.2, and suppose (n — A) -1 <oc< (n— 1) . We define a 
c €°° smoothing function (j)(x) such that (j)(x) = 1 for x £ [—1/2, 1/2], §(x) = for x £ [—2,2] 
and §(x) + (j)(x _1 ) = 1 for all x G R, so 

U(t,x) = - lim y U(kt a ) + ty( — ) i Sm , 7tY e(tk"+xk) = U*(t,x) + limU K (t,x) 
^ K ^°°\jj^ K { \kt a J ) k K^oo 

(5.1) 
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where 

U*(t,x) = - E <|> (kt a ) & -^e(tk n +xk) , (5.2) 

%, k 

k— — oo 

U K (t,x) = - £> (-L) S ^^e(tk"+xk). (5.3) 

K \k\<K \ kt / fe 

We can write Ux{t,x) as the difference of the two expressions 

and will consider only positive values of k; the contribution from negative values can be 
estimated by the same argument. By summation by parts we obtain 



h I *(&)l d \ E e(tk" + (x±y)k) 

Z7U V Af / ^ l/2t a <k<K 



K 

2ni J a\ U 



1 '' '' fl ^f^)) I «(/*" + (x±y)*)d§. 

VV / J l/2f«<A:<E 



1/2f a ./2t*<k<% 



We can bound the remaining sums over k as 

£ e(rt» + {x±i)k)dZ, « {r 1 + <7 _1 + q$- n } 1/N Z, 1+z q e 

l/2t a <k<£, 

by (1.20), where N = 2" _1 . If ^ > m?~ a then we choose g = qj as in the definition of 
& m , a , and the sum is 0(£, l ~ A ' N+e ); similarly for the first term with K in place of £,. On 
the other hand, if £, is smaller than this we take the value of q which the lemma gives for 
M = (m + l)t~ a , and the sum is O (t~ a+aA / N ~ e y On taking the limit as K — > oo we now 
have 

U(t,x) = U*{t,x) + 0(t aA / N - £ ). (5.4) 

6. Bounding the contribution from large frequencies at rational times 

We need to provide an argument for rational times similar to that proved in the previous 
paragraph for irrational times; specifically we will consider t = u/q — y + . The analysis is 
a little trickier, however, so we begin by considering the related sum given by 



n k=-°° \i / K 



where (j) is as used in (5.1). We begin by showing 

U(t,x)=V*(t,x) + 0(q- l/2 ). (6.1) 
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which allows us to consider the finite sum V* in place of U. 
By Fourier inversion 

oo 

VM-ijm />(„) £ *£**.(* + (,+ *)*)*, (6.2) 

— oo 1*1—^ 

where $ denotes the Fourier transform <j>, which can be bounded using repeated integration 
by parts by 

oo 

$ft) = y<K$)e(-$n)d$«n- v ,v = 0,1,2,... • (63) 

In order to change the orders of limit and integral in (6.2) we note that for any y 

l 
1 v-, sin27ty& , „ N f , N sin7i(2A"+l)M 

- Y —L- e (tk n +yk)= U(t,y + u) ^ —du (6.4) 

7i ,,rrf„ k J smnu 

\k\<K 

where U(t,y + u) is piecewise constant, as discussed above, and consider the integral over 
[0, 1.2]; the other half can be treated similarly. The function g(u) = sin(7tw)/7tw is analytic 
and bounded on [0,°°), and its antiderivative G(u) is also bounded on [0,°°) as can be seen 
by observing that it tends to a finite limit as u — )■ oo. Breaking up (6.4) into the integrals over 
subintervals where U(t,y + u) is constant, if [a,b] is any such subinterval then it contributes 

b 

(2K+1) fu(t,y + u)^—g((2K+l)u)du. 
J smiiu 

a 

Integrating by parts the integral in (6.4) is 0{q), uniformly in K and y. The Lebesgue 
Dominated Convergence Theorem can then be applied in (6.2) to change the orders of limit 
and integral, and obtain 

oo 

rM^/twi^^^).)*,. (6.5) 

Applying Theorem 1.4 to the sum over k it can be evaluated as a value of U(t,x), and 

oo 

V*(t,x) = - £ £ e q (u^-vw)f^)u(o,x+- + ^\dt\. (6.6) 

" w(mod<j) v(mod</) _ OQ ^ " " ' 

Using the bound (6.3) the integral over r\ is 

oo 

u(o,x+-j f$(i\)di\+ f Hn)\u(o,x+- + ^J -u(o,x+-j\dr\ + 0{q 



for any A > 0, and applying this in (6.6) the first term gives U(t,x) precisely, by Theorem 
1.4. The integrand in the second term vanishes unless one of the conditions 

-y<x + w/q <y , -y<x + w/q + r\/q 4 < y 
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is true and the other false; this can happen in four ways, each of which implies that w/q 
is in an interval of length 0(q~ 2 ), since the range of integration in r\ is of length q 2 . For 
sufficiently large q there can be at most one w/q in each of these intervals, hence these terms 
contribute O (q~ 1 ^ 2 ) to (6.6). (Here we have used the bound (1.16) for the complete sum 
over v modulo q.) This proves (6.1). 

Recalling that §(kt a ) + §(l/kt a ) = 1, for sufficiently large q we have §(kt a )§{k / q A ) = 
§(kt a ) and the definition of V* and (6.1) imply that 

U M - ,.„,-) + Ij> (jL) * ( * ) ^e ( V + *) + „(,-/») . 

To estimate the second term it is sufficient to bound the sum 

where y = x± y. Summing by parts in a similar fashion to section 4, 

1 4 

but rather than apply (1.20) as we did for irrational t, we use (1.19) and obtain 

2q* 
W(t,x) « J ft' 1 /" + ^-n/N^n/N-l/N q -l/N +l) @/N + q W)logq} £,-*<% 

\/2t a 
^[t a l N + q- l l N + t m l N q l l N ]q\ 

Since oc < \/{n — 1) < 1, we have t a > t > q~ } so the second term can be dropped, as 
can the error term from (6.1). Furthermore the second term dominates the third only if 
t < q-^-M"- 1 ) < q- 1 , which is false for sufficiently large q, so t a > jVC"- 1 ) > q- 1 /^ 1-1 ), 
and 

U{t,x) =U*{t,x) + 0(t an / N qV N+£ ) . (6.8) 



Supposing now that t < mq 1 / (X ( n A ) this error term is 0(t aA l N e ) as in (5.4) and Theorem 
1.5. 

7. An asymptotic estimate for U*(t,x) 

To estimate U* we write it as 

x+y 



U*(t,x)= [ ]T §{kt a )e{tk n + r[k)dT[ 



x-y k ~ 
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and consider the inner sum. Note that for this part of the analysis it makes no difference 
whether t is rational or irrational. Applying the Poisson summation formula 

CO oo 

^(kne(tk n + n k)=^ U(yne(tf + (l\-m)y)dy = t- a Y j U{y)e{F m {y))dy 

k rn m ** 

— CO — oo 

where F m (y) = t l ~ na y" + (r\- m)r a y. For m ^ 

|T|-m| > |m|-(|x|+y) > l-(|x|+y) >0 
so for sufficiently small t 

|^(y)|>|m-Ti|r a -itf 1 - w br 1 >^|m-T,|r a >i(H-(W+Y))r a 

since a < l/(n — 1). Thus for m/0, F^(y) does not vanish in the region of integration and 
we can integrate by parts twice and find 



co oo 

f**w» = -&f£{Tk£{M>}h™»« 



t 2a 



(H-(W+y)) 2 



so the contribution from all m / is thus 

x+y 



r a f £ f^y)e(F m (y))dy^t a ^(\m\-(\x\+y)r 2 ^t a , (7.1) 



x -y m^0_ 



where the implied constant depends on y but not on x, and hence 



x+y 
-all ^(,.\„(A-na,.n 



U*(t,x)=r a J J §{y)e{t v - na y n + T\r a y)dydT\ + 0{t a ) 
x— y — oo 
(x+y)t- l ' n oo 
= J j \{yt a - l l n )e{f+T\y)dydT\ + 0{t a ). (7.2) 

If wis even, the equation ny"~ x +r\ = 0has a unique solution yo = (—T\/n) 1 '( n ~ 1 '; if nis 
odd there at most two solutions, given by ±yo if present. In either case the solutions are in 
the range of support of §(yt a ~ l l n ) if and only if |r|| < w 2"~ V~ 1 /"~ a( "~ 1 ). If x > ythen for 
sufficiently small t we have (x — y)t~ l l n > n 2 n t x ~ 1 l n ~ a ^ n ~ l \ since a < \/(n— 1), so there 
are no solutions yo hi the support of §(yt a ~ l ' n ), and 

|n/- ! +r|| > ri - ny"' 1 > ^-,an-i t i-i/n-a{n-i) >x\- ^p-r 1/n > x\/2. 
Integrating by parts twice in y we can now bound the integral in (7.2) by 

-JL / M (_■_ ' («ggg> Ut^ «,- (j -y)-' . 
47i 2 J J dy 1 ny n - l +r\ dy | m^+Ti | | v ' ; ' v " 

|-„ »A.-l/n-oo ^ ^ ) ) 



(x-y)t 
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A similar calculation shows that if x < — y then the integral is 0{t a \y+x\~ l ); we can com- 
bine these two cases with (7.1) to give U*(t,x) <C t a \x 2 — y 2 ] -1 where the implied constant 
depends on y 

On the other hand, if — y < x < y then this same reasoning applies to the two tails 
r) > (x + y)^ 1 /" and r| < (x-y)^ 1 /" so that in this case 

co oo 

U*(t,x) = J f Wyt a - 1 ' n )e(?> + Tv)dydT\ + 0(t a \x l -'f\- 1 ) . 



— CO — CO 



The inner integral in the remaining expression is the Fourier transform of ^(yt a ~ 1 ^ n )e(y n ), 
hence the outer integral is inverting the transform and gives the value of the original function 
at 0, so combining these estimates with (7.1) we have 

n*(tr) = i l + ^ x2 - y1 ^ if|x|<Y 
{,) \0{t*\x 2 -f\-') if|x|>y 

and combining this last with (5.4) and (6.8) we obtain the first claim in Theorem 1.5. 

The transition between these two cases is the main concern in Theorem 1.5. We will 
consider x near y omitting the details for x near — y which are very similar. If we renor- 
malise values of x as x = y+st x l n for s in some fixed interval [—5,5], then (7.2) becomes 

s+2f- l l" oo 

U* (t,y+st 1 ^ = J J\(yt a - l l n )e{y n +T\y)dydi\ + 0(t a ). 

S — oo 

Since \/n < a < l/(n — 1) by hypothesis, — 1 < (1 — om)(n — 1) < 0, and hence for r\ > 

s + 2yt- l l n 

|„/- 1 +T 1 |>r 1/ "-? (1 - (W)(n - 1)/ ' ! »r 1/n . 

Integrating by parts twice as above, 

co oo 

J J t?{yt a - l ' n )e(f + x\y)dydx\ < t a 
and hence 

co co 

U*(t,y+st l l") = f f §(yt a - l l n )e(y n +i\y)dydi\ + 0(t a ) . 

g — oo 

Suppose n is even. Were the integral in r\ over (— oo^oo) then the double integral would 
be the Fourier inverse at zero of the Fourier transform of §(yt a ~ l l n )e(y n ), so the double 
integral would be equal to 1. Since the integral over r\ € (0,°°) is half this, 

U* {t 1 y+st l ' n \ =\~ f U{yt a - l/n )e{y n +r [ y)dydr [ + 0{t a ) 



\~~k,\ *c>* a " 1/, w) &in2ns y - +0 ^ 
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The calculation is completed by noting that by integration by parts 



' () 

and similarly for y" — sy, and hence 



y J dy I y (ny n 1 + s) 

n v 



U* (t,y+st l l n \ = f e(y n )sm2nsy^ + O (t an - 1 ) . 



On the other hand, if n is odd then (7.2) becomes 

OO CO 

U^t^+st 1 ^ = f U{yt a -V n )cos2n{y n +r}y)dydr\ + 0{t a ) 

S — oo 

S °° 

l -- J J \(yt a - 1 /")cos2n(y n + r\y)dydr\ + 0(t a ) 



2 

o -°° 

CO 

= --— /'^(^ a " 1/ ")sin27l(/+ri3;)^ + 0(? a ) 
2 271 J y 

— oo 

which gives the form in Theorem 1.5. 
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